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Relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts 
are widely considered to be magnetohydrodynamically driven by black hole accretion 
systems, although conversion mechanism from Poynting into particle kinetic energy flux 
is still open. Recent detailed numerical and analytical studies of global structures of 
steady, axisymmetric magnetohydrodynamic (MHD) flows with specific boundary con- 
ditions have not reproduced as rapid an energy conversion as required by observations. 
In order to find more suitable boundary conditions, we focus on the flow along a poloidal 
magnetic field line just inside the external boundary, without treating transfield force 
balance in detail. We find some examples of the poloidal field structure and correspond- 
ing external pressure profile for an efficient and rapid energy conversion as required 
by observations, and that the rapid acceleration requires a rapid decrease of the exter- 
nal pressure above the accretion disk. We also clarify the differences between the fast 
magnetosonic point of the MHD fiow and the sonic point of de Laval nozzle. 

Subject Index xxxx, xxx 



1. Introduction 

Collimated outflows, or jets, with relativistic speeds are observed in active galactic nuclei 
(AGNs) and galactic microquasars, and are presumably driven in gamma-ray bursts (GRBs). 
A widely discussed model for these is that rotational energies of an accretion flow around 
a black hole (BH) and/or of a BH itself are magnetohydrodynamically converted into an 
outflow energy [1, 2], in analogy to pulsar winds [3]. (Alternative possibility that thermal 
energy of a BH accretion flow is transferred into an outflow is also discussed [e.g., 4-6].) In 
relativistic jets driven magnetohydrodynamically, Poynting flux dominates the total energy 
output at the launching site, while observations of AGN jets strongly suggest that the energy 
is dominated by particle kinetic energy at the emission site [e.g., 7-9]. It has been actively 
debated how Poynting flux can be efficiently converted into kinetic energy flux. 

In oder to discuss the energy conversion, ideal magnetohydrodynamics (MHD) is the sim- 
plest approximation. We consider steady, axisymmetric outflows as further approximation.^ 
In this context, the poloidal velocity of the flow has to exceed fast magnetosonic speed (i.e., 
pass the fast magnetosonic point) to continue accelerating towards infinity [13]. At the fast 
magnetosonic point (or the fast point), the bulk Lorentz factor T becomes ~ 

for £• > 1 



^ Unsteadiness and/or non- ideal MHD effects can be important for energy conversion from Poynting 
into kinetic energy flux [e.g., 10-12]. 
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where E is the total energy flux per unit rest energy flux (more precisely, see Eq. 31 below) 
[13-15]. Thus a major energy conversion from Poynting into kinetic energy flux, up to F f , 
should occur beyond the fast point. Self-similar analysis of full MHD equations including 
the transfleld force balance has shown that an efflcient acceleration beyond the fast point 
can occur but that a major energy conversion requires very long distance [16, 17]. 

Pioneering works by Komissarov et al. (2007, 2009) [18, 19] have performed numerical 
simulations of time-dependent, axisymmetric MHD equations over much extended region, 
obtaining steady solutions of relativistic jets with specific boundary conditions [see also 20]. 
They show that the outflows which are confined by a paraboloidal wall accelerate efficiently. 
In those outflows, poloidal magnetic fields near the jet axis are self-collimated, leading to 
a decrease of the field strength and an efficient acceleration in the main body of the flow. 
This result is supported by some analytical treatments [21-23]. In particular, Lyubarsky 
(2009, 2010) [21, 22] constructed asymptotic solutions of the outflow structure far beyond 
the light cylinder, confined by external pressure decaying as a power law of height, under 
the assumption that Poynting flux is not signiflcantly converted into kinetic energy flux near 
the fast point. 

These studies show a rapid energy conversion up to the equipartition level, i.e., the ratio 
of Poynting to kinetic energy flux a ~ 1, but after that, the conversion rate dramatically 
decreases. Then the distance where cr < 0.1 is realized is too large, compared with the obser- 
vational suggestions that a < 0.1 at the distance of the emission site ~ 10^ — 10^ rg, where 
rg is the Schwarzschild radius of a central BH [7-9] . 

However, the previous studies [18-22] have considered only special types of boundary 
conditions. Thus it may be useful to clarify whether there can be boundary conditions and 
magnetic fleld configurations leading to more efficient acceleration that is compatible to the 
observational suggestions. Fendt h Ouyed (2004) [24] have raised a hint for this problem by 
only solving the poloidal direction of motion without treating the transfield force balance 
[see also 25]. They show that a rapid energy conversion up to F f near the fast point 
occurs if the poloidal field strength evolves along a field line as a power law of cylindrical 
radius like Bp oc r"^""? with g > 0, although the corresponding structure of the magnetic 
field was not specified. In this paper, we investigate magnetic field structure required for a 
rapid major conversion of energy. We focus on the fluid motion along a poloidal fleld line 
just inside the external boundary (see Figure 4 below) and show some examples of the field 
line shape and the external pressure profile for a rapid energy conversion (Section 4) . We do 
not treat the transfield force balance in detail. Before showing our findings, we review basic 
theory of steady, axisymmetric, relativistic MHD fiow (Section 2) and discuss properties of 
the fast point a bit more deeply than previous papers (Section 3). 

2. Basic Equations 

We make a brief review of theory of steady, axisymmetric, relativistic ideal-MHD flow [13, 
17, 21, 26-28]. We also assume that the flow is cold, and neglect effects of gravitational field, 
for simplicity. Basic equations are then equation of motion 



pc^(u • V)u = [(V • E)E + (V X B) X B] , 



(1) 



equation of continuity 



V • (pu) = 0, 



(2) 
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Maxwell equations 

V • B = 0, (3) 

V X E = 0, (4) 

and ideal MHD condition 

E+-vxB = 0, (5) 
c 

where u = Fv/c, and T is Lorentz factor of the flow. 

In the axisymmetric configuration, it is convenient to divide the magnetic field into poloidal 
and toroidal components, B = Bp + i?<^e<^, where the poloidal component can be written as 

Bp = ^V*(r,z)xe^. (6) 

We have adopted the cylindrical coordinate system (r, (p, z). Prom Eq. (6) we have B • V^' = 
0, i.e., ^' is constant along each magnetic field line, ^(r, z) is called flux function, representing 
the toroidal component of the vector potential times r or the total magnetic flux penetrating 
within a circle of r for a given z. 

Eq. (4) and Eq. (5) guarantee that E^^ = and E ■ B = 0, respectively. Then the electric 
field can be written as 

E = -^nV*. (7) 

The toroidal component of Eq. (4) leads to B • VO = 0, i.e., J7 = is constant along 

each magnetic field line. We also have a form E = (rOBp/c) x e^. 

Dividing the particle velocity field into poloidal and toroidal components, v = Vp + v^e^, 
the ideal MHD condition gives us 

Vp = KBp, (8) 
v^- KB^ = rn{^), (9) 

where k is a function of (r, z) . 
Eq. (2) together with Eqs. (3) and (8) reduce to 

AnpTK = r/(*), (10) 

where rj = ri{^) is constant along each magnetic field line and represents the mass fiux per 
unit magnetic flux. 

Finally we divide Eq. (1) into the direction parallel to Bp, the azimuthal direction, and 
the direction normal to Bp. The component parallel to Bp is 

d sin^ 2 -^v ^ / N /-■-■% 
^^9/^^ = ^"--4^9Z^'^^)' 

where I is the coordinate of the direction of Bp, and 6 is the angle between the / and z 
directions, i.e., sin^ = dr/dl. The azimuthal component is 

which reduces to 

where J0.{^) is constant along each magnetic field line and represents the total angular 
momentum flux per unit mass flux. Differentiating Eq. (13) with I and substituting it into 
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Eq. (11) together with Eq. (9), we obtain 
which means 

-^-<*). (-) 

where <f (^) is constant along each magnetic field line and represents the total energy flux 
per unit rest energy flux. The ratio of Poynting to kinetic energy flux is given by a = 
—rBipQ./{(?rjr) = {£ — V)/T. The component of the equation of motion normal to Bp, i.e., 
transfield force balance equation, is given in [21], which will not be treated in detail in this 
paper. 

2.1. Bernoulli equation 

If Bp{^,r) is given, one can solve p,T,K.,v^, and B^ as functions of {^,r) with 
ri{^),Q{^),C{^), and as parameters. Eqs. (9), (10), (13), (15), and a relation 

— Up — u"^ = 1 reduce to an equation for Up as a function of x = r/ric for a fixed 

^2 , 1 _ c2 ^^(l - 4 - - [xHl - 4) - M^xl]' 

P ^2(1 - - M2)2 ' ^ ^ 

which is called Bernoulli equation [13]. Here ric = c/Q is the radius of the light cylinder. 



x. = ^|p (17) 

is the radius of the Alfven point divided by ric, and 

M = (18) 

is the Alfven Mach number. Derivative of the Bernoulli equation with x has a critical point 
where the poloidal 4-velocity equals to fast magnetosonic speed, 

2 B^pjl - x') + B^^ _ , 
= ^2 = «/■ (19) 

This is called the fast point [13]. A solution of Up growing towards infinity has to pass both 
Alfven and fast points. We call this a wind solution. Since we have assumed that the flow 
is cold, slow magnetosonic point does not appear. Eq. (16) can be rewritten as the fourth 
order algebraic equation for Up, which can be solved for given Bp{^,x)/r]{'*if), XAi'^), and 
£{^) [13]. Corresponding v^/c and B^/{cr]£) can be calculated by 

^ x\l-x\)-M^x\ 
c x{l-x\-M^) ' ^ ^ 

_ X — Xj^ 

cqS x(l-x2-M2)' ^ 

For given Bp{r), rj, Q, and C along a field line, a unique value of £ is determined for a 
wind solution. More physically, when Bp(r), ri, 0, and the velocity at the inlet Up^in are given 
along a field line, Eqs. (9), (13), (15), and — Up — u"^ = 1 evaluated at the inlet, together 
with the condition of the fast point, determine the values of £, and £ and the solution of 
Up{r). 
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If the solutions for all the range of ^ satisfy the transfield force balance equation and 
certain boundary conditions, they are a solution of the whole outflow structure, although 
we will not treat the whole structure in this paper. 



2.2. Asymptotic relations 

Asymptotic forms of various quantities beyond the light cylinder to the lowest order are 
easily found along each magnetic field line [17, 21]. Eliminating from Eqs. (13) and (15), 
we obtain 



c rn \ r ^ 

Wc may assume that Vy,,in <C c and Hn^^ < c at the inlet, and then we obtain £ — ~ 
Tin. Thus we may write 

!V = ^ _ £i£ y (23) 

c rn\ r J ^ ^ 

For a wind solution in which F S> 1 at r S> He, we have oc r~^. 
Since Vp ^ c and <^ c at r ^ ric, we have from Eq. (9) 

^ - ? » 1- 

-Dp C 

In other words, the poloidal velocity dominates the total velocity, while the toroidal field 
dominates the total magnetic field beyond the light cylinder. Eqs. (15) and (24) give us 

£ ^r + ^B^, (25) 
rjc^ 

for r ^ ric- This equation reads that F grows when BpV^ decreases along a field line beyond 
the light cylinder. This can be also understood by Eq. (11), in which the last term of the 
right-hand side represents the Lorentz force. This force accelerates the fluids when |-Bi^|r 
(oc Bpr^ at r S> He) decreases along a field line. 

For monopole magnetic field structure, i.e., Bpr^ = const, the flow cannot be efficiently 
accelerated. More specifically, in this case, the fast point is located at infinity [13, 14]. 

For later calculations, we show two useful equations for r S> ric here. Eliminating Vp and 
from the relation {vp/cf' + {y^jcf' + l/F^ = 1 by using Eqs. (9) and (23), we obtain 



Bl _ 
Bl ~ c2 



1 c2 „ / 1 c2 



l + -,-^ + -J, 



F2 r2J72 ^ V F4' r2j^2r2 ' ^4J^4 



(26) 



We also have 



, \ C2 2c2 F,, „ / c2 



c 



,2 r2 r2j^2 ^1^2 F V r2n2r2 } 



+ 157^^ + ^? -57^ • (27) 



2.3. Wind solutions 

Fendt & Ouyed (2004) [24] solved Eq. (16) for cases of Bpv"^ oc r^'^ with q > 0. Following 
this formulation, we also solve Eq. (16) and show a wind solution for q = 0.1 in Figure 1. 
The parameters are given as Bp{x = l)/(cr/) = 100, xa = 0.995, and £ = 91.6. Note that if 
Vip <^ rCl and T <^ £ at r = ric (i.e., x = 1), we have £ ~ Bp{x = l)/cr] from Eq. (15). The 
wind solution describes the flow passing through the Alfven point at x = xa and the fast 
point at X ~ 3.2 and accelerating towards infinity. Our result is consistent with that in [24]. 
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Fig. 1 Solution of the Bernoulli equation with BpX^ oc x~'^ with q = 0.1. The parame- 
ters are Bp{x = l)/{cri) = 100, xa = 0.995, and £ = 91.6. The dashed line represents the 
maximum level of the Lorentz factor, T £. 




Fig. 2 Solution for Vip/c (left) and \bip\ = \B^\/{cr]£) (right; solid line) and bp = Bp/{cr}£) 
(right; dashed line) corresponding to Figure 1. 

As can be understood by Eq. (25), the acceleration rate dT/dx monotonically decreases for 
cases of Bpr'^ oc r~^. This behavior is clearly seen in Figure 1. 

Other quantities, v^/c and B^, corresponding to Figure 1 are shown in Figure 2. We can 
see that oc and \B^p\/Bp ~ x at x ^ 1, as discussed in the previous subsection. 

For larger q, the energy conversion is more rapid, as demonstrated in [24] (which we 
also confirmed). Therefore for cases of Bpr"^ oc r~^, we can have an arbitrarily efficient 
acceleration for large q. 

3. Magnetic nozzle 

In this section, we discuss behavior of the flow around the fast point a bit more deeply than 
previous papers. Komissarov et al. (2010) [29] derived an asymptotic equation for r S> He 
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from Eq. (11), which is equivalent to 



^s-r J s 

where 

S = By. (29) 

This equation reads that in order to have dr > 0, 5 has to increase in the region of F < 
{£ - r)V3 and vice versa, and T = {£ — F)^/^ at the fast point. This behavior is quite similar 
to de Laval nozzle for a steady, one-dimensional, non-relativistic hydrodynamic flow [30] , and 
often called magnetic nozzle effect [15], where 1/S corresponds to the cross section of the 
magnetic nozzle. Some results of the numerical calculations by Komissarov et al. indeed show 
that S has a maximum (sec Figure 10 of [19]). On the other hand, in the wind solution shown 
for the case of the monotonic decrease of 5, -Bpr^ oc r^'^ with g > 0, monotonically grows 
from non-relativistic to super-fast speed, passing through the fast point, as first demonstrated 
by [24] and also by Figure 1. This does not appear consistent with Eq. (28). 

In order to solve this apparent inconsistency, we rewrite Eq. (11) into a similar form as Eq. 
(28) but retaining the next order terms. As shown in Appendix, Eq. (11) can be rewritten 
by using Eqs. (9), (14), (26), and (27) as 

dV ^ c^Fin dV f _^ ^ 2Fin^ {—v^pdvip) —dS 



Note that this equation is valid only for r > r\c. Clearly, this equation is not exactly the 
same as the case of de Laval nozzle. Even at and inside the fast point (i.e.. Up < uj), the 
second and third terms of the left-hand side allow dF > while dS < 0. Thus the monotonic 
decrease of S discussed in [24] can have a solution of the flow passing through the fast point. 
For r ^ ric, however, the second and third terms of the left-hand side are negligible, and 
uj ~ B'^/iAnpc^T^) ~ (f - F)/F, and then we have Eq. (28). At the fast point, we have 

F-(^-F)^/3 r = rf. (31) 

This property has been already shown in [15]. 

Now let us examine which shape of S leads to more efficient acceleration, the monotonic 
decrease shape or the nozzle shape. In order to do this, we compare wind solutions of Eq. 
(16) with different models of S. The first model is 5 oc with g = 0.1 as shown in Figure 
1, while the second one is a broken power-law form that has a maximum point, 

with p = 0.5, q = 0.1, s = 3, and C = 2.0 (see Figure 3 left). For these parameters, S has a 
maximum at x ^ 3.4. We set Bp/(cr]) as taking the same value as that of the first model at 
a; = 30, so that Bp is almost the same for the two models far outside the fast point. We set 
the same value of xa = 0.995 for the two models. 

Figure 3 (right) shows the two wind solutions. The first model, represented by the solid 
line (which is the same as Figure 1), has £ = 91.6 while the second model, represented by 
the dashed line, has £ = 89.2, a bit smaller than the first model. The radius of the fast 
point is a bit larger in the second model. However, the distances spent for the major energy 
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Fig. 3 Different structures of S along a magnetic field line (left) and corresponding solu- 
tions for the Bernoulli equation (right). The solid line represents the model of <S oc x~^-^ 
(same as Figure 1) and the dashed line represents that of Eq. (32). 



conversions are very similar. This means that the total acceleration efficiency hardly depends 
on the field structure inside the fast point. In summary, one has F ~ (iS — F)^/'^ at the fast 
point in general, and the total acceleration efficiency is determined by the poloidal field 
structure outside the fast point. 

4. Magnetic field shapes for efiicient acceleration 

In Section 2, we have seen that an arbitrarily efficient acceleration is available in the case of 
Bpr^ oc with q> along a field line. However, two-dimensional structure of the poloidal 
field lines that gives such a scaling of Bp{r) is not clear. In this section, we show some types 
of the flux function ^{r, z) for a rapid acceleration just beyond the fast point. 
Before doing it, we review an instructive example of the flux function that is given by 

y+ai^) = Amx-pi^)r^^\ (33) 

where y = z/r\c and = ^/^o- This form of * has been discussed in the literature [31, 32]. 
We assume (' > 0, A' < 0, p' > 0, and a' < to guarantee that ^' is smaller for more inside 
field lines, where the prime means the derivative with tp. By calculating Br = {—l/r)d^/dz, 
B^ = {l/r)d^/dr, and Bp = y^B^f^Bl, we have 



= ^ " ^ " 




(34) 



If a > 1, the last factor has an effect of decreasing BpX^ along the field line for the range 
of {x — p)/(y + Cj > i.e., for the range of B^. > B^- In the second factor, C /{v + C) and 
ap' /{x — p) increase Bpx'^, A' /A does not change Bpx'^, and — a'ln(a; — p) decreases Bpx"^ . 
Thus only for a' 7^ 0, a long-lasting decrease of Bpx'^ is realized, although it decreases only 
logarithmically. Below we discuss a similar type of the flux function with a' / as Type 1. 

^.1. Model and parameters 

Poloidal magnetic flux within a BH accretion system is not infinite in reality. We assume 
that the system has a finite region with non-zero poloidal magnetic flux around the axis, and 
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Fig. 4 Schematic picture of our model for Type 1 and 2. The whole outflow region with 
non-zero poloidal magnetic flux is confined by the external medium with thermal pressure 
and/or magnetic pressure of closed fields. We fix the shape of the external boundary as the 
solid line, and discuss the flux function for the region just inside the boundary, *o — < 
'^{r,z) < and wind solutions for the field line of ^'o- The dashed curved line represents 
y = Aqx"'° , where y = z/ric and x = r/r\c. 

this region drives a relativistic outflow. The finite magnetic flux is defined as ^'o (see Figure 
4). The flow is confined by the pressure of the external medium, which includes thermal 
pressure and/or magnetic pressure of the closed fields. In this case the pressure balance 
condition should be satisfied at the boundary [cf., 21]. Since the magnetic field measured at 
the plasma rest frame is B' = {B^ — the boundary condition is given by 



We fix the shape of the external boundary, and discuss the flux function for the region just 
inside the external boundary, — 6"^ < ^{r, z) < ^o, and wind solutions for the fleld line 
of ^0- 

We assume that the boundary at r > ric has a shape of 



(i.e., ^'(r, z) = *o is equivalent to Eq.(36)), where and oq are constants (except for the 
final case in Type 2 below). At r < ric the boundary and the last field line are assumed to 
deviate from this shape and be anchored to a thin disk at tq with Keplerian rotation, as 
illustrated by the solid line in Figure 4. However, for simplicity, we use wind solutions of the 
Bernoulli equation along the field line shaping as y = Aqx""" for < a; < oo (i.e., the dashed 
curved line in Figure 4). This treatment is justified since the acceleration efficiency does not 
strongly depend on the field structure inside the fast point, as discussed in Section 3. For 
wind solutions we obtain, Pext{z) can be deduced by Eqs. (35) and (36). 




(35) 



y = Aqx' 



(36) 
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The parameter is given by 




(37) 



Then the hght cyhnder radius is 



nc(*o) 



c 




(38) 



For the field Hne with Ro = 4, as an example, we have r\c — 10 rg. 

The wind solution Up{x) and corresponding external pressure profile Pcxt(2;) depend on 
the parameters Aq and ao, but not on i?o, which just converts (x, y) to (r, z) by using 
Eq. (38). To have rough constraints on Aq and ao, we may use a simple opening angle 
estimate at the emission site as Ve/ze ~ 0.1. Substituting this relation into y = Aqx^° with 
ao = 2, we have X 10^ ~ AQ^Ze/rs). Observations suggest that Ze ~ 10^ — 10^ rg. For 

Ze/fs = 10^, we have ~ 0.4, 1, and 3 for Rq = 2, 4, and 8, respectively. For Ze/vg = 10^, we 
have Aq ~ 0.04, 0.1, and 0.3 for Rq = 2, 4, and 8, respectively. For cq = 1, we have Aq ~ 10 
for any value of Rq. Below we will adopt ylo = 0.1 and ylo = 1 with ao = 2 and for Aq = 10 
with ao = 1 as the fiducial values, for which we show calculation results. 



First let us consider a type for which the flux function for ^'o — (5* < ^' < ^'o is given by 



where a! < and a^ip = 1) = cq. The parameters D and d are constants. This type is sug- 
gested to make a long-lasting acceleration from the discussion below Eqs. (33) and (34). The 
constant d should be ^ 1, since all the magnetic field lines intersect and the field strength 
would diverge at x = d. For the field line of ^' = *0) we have 



where we have defined aQ = a'{ip = 1). 

Figure 5 (left) shows a wind solution for BpX^ given by Eq. (40). The parameters are 
given as Bp{x = l)/{cr}) = 100, d = 10"^, D = 10"^, ao = 2, = 0.991, and S = 55.8. This 
parameter choice corresponds to the field line shape of ^' = ^'q as Eq. (36) with gq = 2 and 
Aq = 0.1. For x/y > oq, i.e., x < 5, Bpx'^ rapidly decreases, roughly scaling as oc x/y oc x~^, 
i.e., BpX^ oc x~'^ with g = 1, so that the acceleration is very rapid. This causes the fast point, 
where F ^ £^/^ ~ 3.8, to be closer to the Alfven point than the case of Figure 1 (g = 0.1). For 
X > 5, however, the last factor of Eq. (40) is roughly constant, so that Bpx"^ decreases only 
logarithmically. As a result, wc have a < 0.2 at x > 170, corresponding to y > 2.9 x 10^ and 
then z > 4 X 10^ RQ^'^r^, although a < 0.1 is realized only at x > 10^, corresponding to z > 
V2 X 10^i?o^Vs, which is very large distance compared with the observational suggestions. 

To illustrate the magnetic field structure of this type, we plot the field line of ip = 0.95 
with the dashed line in Figure 5 (right). This is the line of Eq. (39) with the power law index 
set to be a('0) = ao + ao(V' ~ !)• We have assumed that *o/(c'7''ic) = 10 x In(lO^) 2± 46 and 
given Qq by using Eq. (40) evaluated at x = 1 with Bp{x = l)/{cr)) = 100. 



4.2. Type 1 




(39) 




(40) 
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Fig. 5 Left: Solution of the Bernoulli equation with BpX^ given by Eq. (40) for Up. 
The parameters are Bp{x = l)/{cri) = 100, d = 10"^, D = 10'^ oq = 2, xa = 0.991, and 
6 = 55.8. The dashed line represents the maximum level of the Lorentz factor, T ^ 6. Right: 
Poloidal field lines of ip = 0.95 for the flux function given by Eq. (39) (dashed line) and Eq. 
(42) (dot-dashed line) for which ^Q/{crirf^) 2± 46 is set. The solid line represents the external 
boundary ■0 = 1. 




0.1 1 10 100 1000 10000 



X 

Fig. 6 External pressure as a function of x corresponding to Up in Figure 5 (left) . 



Figure 6 shows the external pressure calculated by Eq. (35). Pcxti^) decreases very steeply 
just beyond the light cylinder, and it roughly scales as oc x"^ far outside the fast point. This 
behavior can be understood by an approximate equation of Eq. (35) , 

P^^'^^^f^' (41) 

which is derived by using Eq. (26) for a; > 1. For the rapid acceleration region, 1 < x < 3, 
the external pressure Pext oc S^r^/F^ oc x~^/T^ decreases by a factor of ~ 10""^. For x S> 5, 
F ~ const, and the external pressure scales as Pext oc Bpv'^ oc x~^/{\nx)^. Since z oc x^, we 
have Pext oc z~^/{\nz)^. 
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Fig. 7 Solution of tiie Bernoulli equation with Eq. (44) for Up (left) and Pexti^) (right). 
The parameters are Bp{x = l)/icr]) = 100, = 0.1, Oq = 2, Fq = 0, 6 = 3, xa = 0.990, and 
£ = 48.2. The dashed line represents the maximum level of the Lorentz factor, F Ri 



4.3. Type 2 

Next we consider a flux function given by 

y = Aqx'"' + F{'iIj)x\ (42) 

for *o - < * < *o, where uq < b, F{il)) < Ao, F' < 0, and b is constant. For this flux 
function we have more efficient and rapid acceleration. Eq. (42) provides 



2 _ ^0 Aqx""-'' + F IfxY , f Apaox^o-^ + Fbx^-^ , 
-^2 (_^,) \l[y) +[ Aox-o-i+Fxb-^ J ■ ^^"^^ 

Choosing the boundary shape as Eq. (36), we set Fo = F{il) = 1) = and obtain for the field 
line of ^'o 

where we have defined F^ = F' {ijj = 1) . 

Figure 7 (left) shows a wind solution for BpX^ given by Eq. (44). The parameters are Bp{x = 
l)/(cr/) = 100, ^0 = 0.1, ao = 2, 6 = 3, rcyi = 0.990, andf^ = 48.2. In this case, BpX^ decreases 
as a power law function of x, not logarithmically (i.e., BpX^ oc x""^ with g = 6 — ao = 1), even 
outside the rapid decrease phase of x/y > ao. As a result, we have cr < 0.1 as early as a; > 6. 
This corresponds to z > 3.6 x Ts ~ 40 rg for Rq = 4, which is small enough to obtain 

0" < 0.1 at the emission site, as suggested by observations. 

The dot-dashed line in Figure 5 (right) represents the poloidal field line oiip = 0.95 for this 
type of the flux function. This line is plotted as Eq. (42) with F{tp) = Fq + Fq{'iIj — 1) and 
with the same value of '^o/{cr]rf^) as that for Type 1, i.e., the dashed line in Figure 5 (right). 
This clearly shows that the field lines of Type 2 expand sideways (and BpX^ decreases) more 
rapidly than those of Type 1. 

Figure 7 (right) shows the corresponding structure of the external pressure. It drops very 
rapidly at 1 < a; < 3 and decays as a power law of a; at a; 3> 5, similarly as Case 1. At x 3> 5, 
Eq. (41) implies that Pgxt oc B^r'^ oc oc z"^. 

We also obtain a wind solution for the case of Aq = 1.0 and show the result in Figure 8. 
For this parameter value, the rapid acceleration phase x/y > qq (i.e., x < 0.5) is so short 
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Fig. 8 Solution of the Bernoulli equation with Eq. (44) for Up (left) and -Pext(^) (right). 
The parameters are Bp{x = l)/(cr/) = 100, Aq = 1.0, uq = 2, Fq = 0, b = 3, xa = 0.991, and 
f = 55.8. The dashed line represents the maximum level of the Lorentz factor, T ^ 6. 
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Fig. 9 Poloidal field line of i/^ = 0.95 for the flux function given by Eq. (42) (dashed line) 
for which Bp{x = l)/(cr?) = 100, ao = 1, b = 2, Aq = 10, Fq = 0, and *o/(c??rfJ 2± 46 are 
set. The solid line represents the external boundary ^ = 1. 

that the acceleration is slower than the case of Figure 7. This means that less-collimated jets 
have a slower fluid acceleration in this type of the flux function. Even though the last term of 
Eq. (44) does not strongly contribute to decreasing BpX^, Up increases rapidly just beyond 
X ~ 1. This is simply a generic property for the cases of BpX^ decreasing as a power law 
function of x, as discussed in Section 2.3. While the acceleration is rapid in the initial phase, 
the acceleration rate, dV/dx, decreases monotonically, as can be expected by Eq. (25). The 
range where cj < 0.1 is a; > 18, which corresponds to z > 3.2 x 10^ x y^^o^^^'s ~ 4 x 10^ Ts 
for Rq = 4. This is comparable to the observational suggestion ^ 10'^ rg. 

Next let us examine the cases for ao = 1, 6 = 2, and Fq = 0. The external boundary in this 
case has the conical shape. The field line of V' = 0.95 for Bp{x = l)/{cr)) = 100, = 10, 
and *o/(c'7''ic) — 46 is shown in Figure 9. The wind solution we obtain for xa = 0.992 and 
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Fig. 10 Solution of the Bernoulli equation with Eq. (43) for Up (left) and P^xtix) (right). 
The parameters are Bp{x = l)/(c??) = 100, Ao = 1.0, ao = 2, Fo = 0.03, b = 3, xa = 0.990, 
and 6 = 48.7. The dashed line represents the maximum level of the Lorentz factor, T ^ 6. 

6 = 58.2 (not explicitly shown) is quite similar to Figure 8. This result is expected, because 
the last factor of Eq. (44) is constant, not contributing to decreasing Bpx'^, and the power 
law index of BpX^ for x, —q = ao — b = —1 is the same as the case of Figure 8. The range 
where a < 0.1, x > 18, corresponds to z > 1.8 x 10 x '-^ 2 X 10^ rg, which is small 

enough to be consistent with the observational suggestions. 

Finally, we consider the case in which Fq / 0. In this case, the boundary does not have the 
shape of Eq. (36) but obeys Eq. (42) with Fq ^ 0. For instance, we adopt the parameters as 
^0 = 0.1, ao = 2, = 0.03, and 6 = 3. We find a wind solution for Bpx'^ given by Eq. (43) 
with Bp{x = l)/{cr]) = 100, xa = 0.990, and £ = 48.7, and the result is shown in Figure 10. 
In this case, decrease of Bpx"^ stops at x ^ Aq/Fq ^ 3, so that the energy conversion saturates 
at (T ~ 0.22. The external pressure far beyond the fast point obeys Pext oc BpV^ oc x~^. 

5. Conclusion and Discussion 

Recently, numerical and analytical studies of global structures of steady, axisymmetric, 
relativistic MHD flows with specific external boundary conditions have been significantly 
developed [18-22]. They all showed that the energy conversion from Poynting into particle 
kinetic energy flux is rapid up to the equipartition level, i.e., until a 1 is attained, but fur- 
ther conversion up to o" ^ 1 is very slow. As a result, the radius where o" < 0.1 is too large, 
compared with observational suggestions that a < 0.1 at the emission site z ~ 10^ - 10"^ rs 
[7-9]. 

In this context, we consider that it is important to clarify whether there can be magnetic 
field structures and boundary conditions for which rapid energy conversion up to cr < 0.1 
occurs. We have fixed the shape of the external boundary of the outflow (see Figure 4) and 
looked for general types of the flux function for the region just inside the boundary leading 
to fj < 0.1 soon after the fast point (without discussing the whole transverse structure of 
the outflow). Although we have solved the equation of motion along the field line shaping 
as z oc r"° for < r < GO (dashed line in Figure 4), those solutions are roughly applicable 
for the more realistic field line anchored at vq < Hc (the solid line in Figure 4), because the 
acceleration efficiency is determined by the field structure outside the fast point, as discussed 
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in Section 3 (see Figure 3). We have found a type of the flux function, given by Eq. (42) 
near the boundary, which describes the magnetic flux tubes expanding sideways rapidly as 
shown by the dot-dashed line in Figure 5 (right) and the dashed line in Figure 9, and leads 
to a rapid energy conversion just beyond the fast point (see Figures 7 and 8). In this case we 
can have a < 0.1 at r < 10^ — 10^ rg, consistent with observational suggestions. For the flux 
function given by Eq. (39), which was often discussed in the literature [31, 32], the flux tubes 
expand slower than the above function, as shown by the dashed line in Figure 5 (right), and 
the acceleration is much slower. 

The poloidal velocity of the fluid is parallel to the poloidal magnetic field line, and the 
fluid is accelerated by the Lorentz force when \B^\r decreases along a field line (see Eq. 11). 
Generally, in the relativistic fluid, Poynting flux is not significantly converted into particle 
kinetic energy flux at the fast point (i.e., F « £^^^), so that the major energy conversion 
has to be realized beyond the fast point. The major conversion up to F f (i.e., a <C 1) is 
realized for the poloidal field structure for which S^r^ continues decreasing along a field line, 
which is accompanied by decreasing |-B<^|r. For the fiux function we found, the magnetic flux 
tubes expand sideways rapidly, scaling as Bpr^ oc with g > 0. In this case an arbitrarily 
efficient acceleration is available for an arbitrarily large q, as clearly demonstrated by [24, 25]. 
Here l/{Bpr'^) causes an effect similar to the cross section of dc Laval nozzle [15, 29], although 
the MHD equation of motion allows the flow to pass through the fast point even for the cases 
of a monotonic decrease of Bpr'^, as we have shown by deriving Eq. (30). 

By using the wind solutions for the field line just inside the external boundary, we have 
calculated the external pressures Pcxt as functions of r (also of z by using the assumed form 
z oc r"") for various cases through the pressure balance condition (Eq. 35). Then we have 
found that Pext has to decay very rapidly just beyond the light cylinder (by a factor of 
~ 10^ — 10^ for r\c <r < 3ric) and asymptote to the power law decay Pext oc B^r^ for all 
the cases that we have calculated, i.e., the cases of the major energy conversion just beyond 
the fast point. Such external pressure profile might be due to thermal atmosphere or corona 
with closed magnetic field loops plus dilute wind above the thin disk, although the origin of 
such dilute winds is not so clear. 

For the cases of BpV^ oc r~'^, the acceleration rate dT/dr monotonically decreases (as found 
from Eq. 25), and therefore a rapid major energy conversion requires the initial rapid accel- 
eration phase, which is accompanied by the rapid decay of Pext- If -Pext maintains a rapid 
decay even far beyond the light cylinder, BpV^ has to decrease with increasing q, and the 
acceleration will be more rapid than the cases of constant q, although the cases of variable 
q are not included in the type of the flux function we have discussed. If Pcxt is constant 
beyond the region of the rapid decay region, the shape of the boundary will be cylindrical 
(i.e., r = const) at the region of the constant Pext, where Bp and F will be also constant. 
That is, the acceleration will stop beyond the region of the rapid pressure decay. 

Komissarov ct al. (2007, 2009) [18, 19] found global solutions above the inlet at r < r\c 
with a fixed paraboloidal wall as the external boundary, by using numerical methods. They 
show that the poloidal fields near the jet axis arc sclf-collimatcd, leading to a decrease 
of BpV^ and an efficient fluid acceleration in the main body of the flow, and call this 
collimation- acceleration mechanism. However, their solutions do not include the cases with 
rapid acceleration or rapid external pressure decay just beyond the light cylinder. Our results 



15/18 



indicate tiiat tlie rapid acceleration requires the rapid external pressure decay, which might 
arise in the numerical calculations if the rigid wall suddenly expands near the light cylinder 
as shown in Figure 4. Such cases can be seen in numerical calculations of the global structure 
by Tchekhovskoy et al. (2010) and Komissarov et al. (2010) [29, 33], in which the opening 
angle of the wall is set to suddenly become larger at some large distance. The numerical 
results show that the poloidal field line expands at the transition point and the rarefaction 
wave propagates inward, leading to additional rapid acceleration. Note that cr <C 1 is real- 
ized near the wall. The authors consider the context of GRB jet breakouts from dense stellar 
envelope region into dilute interstellar medium, and called this rarefaction- acceleration mech- 
anism [see also 34]. Thus our model may correspond to the rarefaction acceleration at the 
jet launching site. 

A highly non-uniform distribution of the poloidal field strength at the inlet might cause an 
additional effect for the energy conversion. For instance, if the field strength is much weaker 
near the jet axis, i.e., around the central BH, the flux tubes could expand towards the jet 
axis. 

The field structure inside the light cylinder docs not strongly affect the acceleration effi- 
ciency, but may affect on the mass injection rate. If the angle between the field line and 
the accretion flow is small, the centrifugal force can easily blow the particles against the 
gravitational force, and the outflow becomes non-relativistic [35, 36]. We speculate that the 
shape of the field line nearly vertical above the accretion flow, as illustrated in Figure 4, 
is effective to suppress mass injection from the accretion flow and allow the outflow to be 
relativistic. Small amount of mass can be injected through electron-positron pair creation 
by collisions of two photons emitted from the the corona [37, and references therein] and/or 
through electron-proton creation by decays of neutrons escaped from the corona [6] . 

The rapid external pressure decay just beyond the light cylinder may arise more easily for 
geometrically thin accretion flows. For geometrically thick accretion flows, such as advection- 
dominated accretion flows [e.g., 38, 39], the outflow region will be confined by the dense 
accretion flow (probably with disk wind), so that the rarefaction acceleration seems difficult. 
The central region of GRB jets within collapsing stars are expected to have high pressure 
[e.g., 40]. Thus the rarefaction acceleration may not be applicable for them, although such 
acceleration is effective at the jet breakout sites from the dense region into dilute interstellar 
medium, as discussed above. For compact stars mergers, the density can be very low above 
accretion flow if a BH is formed [e.g., 41], so that the rarefaction acceleration at the jet 
launching site could be effective [34] . 
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A. Derivation of Eq. (30) 

We derive Eq. (30) by taking an asymptotic limit of Eq. (11), which can be rewritten as 

Updup = -^dr - - — "^dirB^p). (Al) 
r AiToc^r 
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Here all the derivatives represent those in the direction of Bp. 
First, we have from Eq. (14) 



and by differentiating Eq. (9) with dO. = Q 



dr 



1 



—^dvp - -^d{rB^) + ^ dS. 



(A2) 



(A3) 



Combining these two equations with Eq. (Al), we obtain 
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where we have used the relation Vpdvp + v^dv^p = c^dT/T^. 

The parenthesis in the first term can be rewritten by using £ — F = —rQ,B^Bp/ {AirpUpC^) 
as 
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where we have used Eq. (23). By making in the parenthesis, the whole equation reduces 
to 
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This equation can be approximated into Eq. (30) by using Eqs. (26) and (27) and neglecting 
terms of orders of c^/r^ri^F^, c^/r^Cl^, and 1/F^ in the parenthesis { }. 
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